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^ . Abstract. An overview of the searches for gravitational waves from radio pulsars with LIGO 

' and GEO is given. We give a brief description of the algorithm used in these targeted searches 

and provide end-to-end validation of the technique through hardware injections. We report on 
some aspects of the recent S3/S4 LIGO and GEO search for signals from several pulsars. The 
gaussianity of narrow frequency bands of S3/S4 LIGO data, where pulsar signals are expected, 
is assessed with Kolmogorov-Smirnov tests. Preliminary results from the S3 run with a network 
QO ' 01 fo ur detectors are given for pulsar J1939+2134. 
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1. Introduction 

Radio pulsars are interesting sources for currently operating gravitational wave detectors P] 12] ■ 
The positions and spin evolutions of these sources are well known from radio observations which 
reduce the parameter space considerably compared to all-sky searches for unknown neutron 
(3JT)! stars. In addition to simplifying the data analysis problem, the smaller parameter space in 

these targeted searches improves the sensitivity to gravitational waves by lowering the detection 
threshold. 

LIGO and GEO have so far conducted two searches for gravitational waves from radio pulsars 
using data from the SI science run and data from the S2 run [3]. Neither of these searches 
provided a detection but upper limits were set on gravitational wave emission for selected pulsars. 
Both of these analyses searched for gravitational waves at twice the rotational frequency of 
pulsars, where we would expect a signal from an asymmetric pulsar. 

In the SI analysis, two techniques were used to set upper limits on gravitational wave emission 
from pulsar J1939+2134 (the fastest rotating known millisecond pulsar) [3]: a Bayesian time- 
domain method and a classical frequency-domain method. The main result from this SI analysis 
was an upper limit on signals from pulsar J1939+2134 of ho < 1.4 x 10~ 22 with 95% confidence. 

Here we will focus on the time-domain Bayesian method which is more suited for targeted 
sources The frequency-domain statistical technique used in the SI paper is currently being 
used for broad all-sky searches for unknown sources, including Einstein@home [SJE]. 

For the LIGO S2 run, the time-domain search was expanded to include all well-known 
isolated pulsars with putative gravitational wave frequencies above 40 Hz |3|. Using the S2 data, 
multi-detector upper limits were set on gravitational wave emission from 28 pulsars including 
J1939+2134 and the Crab pulsar. The tightest limit on gravitational wave strain came from 
pulsar J1910-5959D with a 95% upper limit of ho < 1.7 X 1(T 24 . At the time of S2 this was the 
lowest upper limit ever set on a pulsar by a gravitational wave detector. The four best limits 
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Figure 1. Best sensitivities of the three LIGO and the GEO 600 interferometers during the S3 
data run. 

on equatorial ellipticity for S2 came from the four closest pulsars J0030+0451, J1024— 0719, 
J1744-1134, and J2124-3358 with 95% upper limits on the ellipticity of 4.8 x 10~ 6 , 8.6 x 10~ 6 , 
8.3 x 10~ 6 , and 4.5 x 10 -6 , respectively. 

The analysis of S3/S4 LIGO and GEO data is currently underway with more sensitive data 
and a larger selection of pulsars than in S2. The main change in our search since S2 has been 
the addition of pulsars in binary systems. The timing for all sources was supplied by the Pulsar 
Group at Jodrell Bank Observatory. 

We are currently analyzing S3/S4 data for gravitational waves from all known radio pulsars 
with /gw > 40 Hz for which the source parameters are known sufficiently well. This breaks down 
to a total of 93 pulsars of which 60 are binaries and 33 are isolated. The improved sensitivity 
of the detectors in S3/S4 promises to give interesting results for several sources. For the Crab 
pulsar, we should be within a factor of a few of the spin-down based upper limit [5]. Figure ^ 
shows the sensitivity of the GEO and LIGO interferometers during the S3 run. 

In this paper we will report on some aspects of the new S3/S4 search for signals from radio 
pulsars. In Section |2] we will give a brief review of the technique used in previous and current 
analyses. We will provide an example of a hardware injection of a pulsar signal in the LIGO 
interferometers during S3 in Section |3J In Section ^ we will show the results of statistical tests 
investigating the gaussianity of the frequency bands of the LIGO S3/S4 data containing the 
radio pulsar signals. As an example, we will show preliminary results in Section [5] for a multi- 
detector analysis for pulsar J1939+2134 using LIGO and GEO S3 data. Brief conclusions will 
be given in Section El 

2. Bayesian time-domain method 

The overall pipeline for the analysis of the S3/S4 data for signals from known pulsars is very 
similar to the techniques used in SI and S2. A brief description of the method is provided here 
with more details available in IE] ■ 

The data from the interferometer are heterodyned with the known phase evolution of each 



pulsar and down-sampled to one data point, B k , each kth. minute. These data points B k are 
essentially a 1/60 Hz band- limited time series centered on the instantaneous frequency of the 
gravitational wave signal at the detector. 

We assume that the noise floor for each band was stationary and gaussian over periods of 30 
minutes. As will be shown later in Section^ this is a reasonable assumption for a majority of 
the pulsars. 

We take a Bayesian approach for the statistical analysis. We use Bayes' Theorem to calculate 
the posterior probability, p(a\{Bk}), of a set of parameters a given the binned data, {B k }. For 
this problem, Bayes' Theorem states that 

P(al{Bk}) = P({B k }) ' (1) 

where a represents the set of parameters that could produce the set of data, {B k }, with likelihood 
p({Bk}\a). Note that a comprises of four unknown parameters: the gravitational wave amplitude 
ho, the polarization angle ifr, the inclination of the pulsar's spin axis with respect to the line-of- 
sight l, and the initial phase of the gravitational wave signal (fro- Our prior beliefs in our set of 
parameters are reflected in the prior probability term, p(a). We use the least informative priors 
for most of the parameters in their respective ranges: (fro uniform over [0, 2n], ip uniform over 
[— 7r/4, 7r/4], and i uniform in cost over [—1, 1], corresponding to a uniform prior per unit solid 
angle of pulsar orientation. 

It can be shown that the likelihood of each 30 minute segment of data, labeled by j, is given 

by 



fc 2(j) 




p({B k }j\a) oc I V \B k -y k \ 2 I , (2) 



where y k is the signal model, B k is the are the data points, and rrij = 30 is the number of 
data points in the segment. Equation |21 is equivalent to a Student's t-distribution with 2m j — 1 
degrees of freedom. We note that the lengths of the stationary segments, mj, can be adjusted 
depending on the performance of the detectors. 

The joint likelihood of all the M stretches of data, taken as independent, is therefore 

M 

p({5j|a)a]>({^b|a)- (3) 

j 

3. Hardware injections 

We use the term hardware injections for simulated signals that make the instrument behave in 
the same way as if a gravitational wave signal was present. During the S3 and S4 runs, several 
hardware injections were carried out that mimicked periodic signals from pulsars in the LIGO 
interferometers. These injections provide end-to-end validations of the search codes and the data 
acquisition pipelines. The phase of the signals is properly modulated to simulate the Doppler 
shift from specific directions in the sky. The signals are also amplitude modulated to reflect the 
antenna pattern of the interferometers. 

These pulsar hardware injections increase our confidence that the timing between sites is 
consistent. This is especially important for coherent multi-detector searches where a large lag 
in timing between the detectors could severely reduce the sensitivity of the search. 

Figure 121 shows the marginalised probability distribution functions for ho and (fro for one of 
the S3 pulsar injections. The recovered parameters for this example agree well with the injected 
signal parameters. Work is currently in progress to make sure all the S3/S4 hardware injections 
are correctly extracted from the LIGO data. There are still some unclear issues with the phase 
of the S4 injections but we do not expect any major problems. 
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Figure 2. Recovered parameters for one of the S3 pulsar hardware injections. The dashed line 
represents the signal parameter injected. 



4. Characterizing the S3/S4 LIGO data 

In this section we report on some statistical tests that were performed to assess the gaussianity 
of the S3/S4 data near the frequencies where we expect pulsar signals. For each pulsar, we 
examined the data in a 1/60 Hz frequency band around the instantaneous gravitational wave 
frequency of the signal at the detector. These data are the heterodyned data set mentioned 
earlier. 

We have applied Kolmogorov-Smirnov (Lilliefors) tests to assess whether contiguous segments 
of 30 minutes of this data are normally distributed for each pulsar. This consists of comparing 
the experimental cumulative distribution function (cdf) for each of the 30 minute data segments 
with a gaussian cdf. The K-S test statistic is a simple method to quantitatively measure the 
difference of these two distributions. 

For each pulsar, the K-S statistic was calculated for every 30 minute segment of data and 
compared to the Lilliefors critical values |Hj ■ The critical value for a sample size of 60 (since the 
BkS are complex) for a level of significance of 0.05 is D = 0.114. By definition, at a level of 
significance of 0.05 we would expect approximately 5% of the data to be rejected if the data was 
indeed normal. Figure [3 shows the distribution of the fraction of 30 minute segments of data 
rejected for each interferometer. For the majority of the pulsars between 5-10% of the data were 
rejected. This is quite reasonable and suggests that the data is relatively well behaved for most 
pulsars. For a few pulsars the S4 data is significantly better modeled by a gaussian distribution 
than the S3 data. This is due to an improvement in the performance of the detectors in S4. It 
has been shown that rejecting the non-gaussian S3 data does not significantly effect the final 
results [Hj. This is not surprising since the non-gaussian data will have a small likelihood. 
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Figure 3. Histogram of the percentage of data which is rejected at a 0.05 significance level 
with a Kolmogorov-Smirnov (Lilliefors) test for normality for each of the 93 pulsars with the S4 
(left) and S3 (right) data. 



5. Preliminary multi-detector analysis between LIGO and GEO 

The full and final results from the S3/S4 LIGO and GEO targeted pulsar search will be presented 
in a future paper by the LIGO Scientific Collaboration and Jodrell Bank Observatory. Here we 
present preliminary results for pulsar J1939+2134 using GEO and LIGO S3 data. Figure |1] 
shows the posterior probability distributions for parameters of a signal from J1939+2134 using 
the S3 data. It is clear that no gravitational wave signals are detected. It can be seen from the 
plot that an 95% upper limit on gravitational wave strain from J1939+2134 using S3 data can 
be set to ~ 5 x 10 -24 . This limit is about a factor of two better than the S2 results jlj and a 
factor of 20 better than the SI result [2]. We can expect further improvement when the S4 data 
is added to the analysis in the near future. 

6. Conclusions 

Upper limits on gravitational wave emission from a selection of 28 isolated pulsars were set 
using the LIGO S2 data jlj. A similar analysis, including binary pulsars, is currently underway 
using the S3/S4 LIGO and GEO data. This analysis will include 93 radio pulsars and is nearly 
complete. 

In this paper we have investigated the gaussianity of the LIGO S3/S4 data in narrow frequency 
bands where the pulsar signals are expected. Generally the data are well described by a gaussian 
distribution as expected. We have presented a preliminary multi-detector upper limit for pulsar 
J1939+2134 using LIGO and GEO S3 data which should be significantly improved when the S4 
data is included. 
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Figure 4. Preliminary posterior probability distributions for pulsar J1939+2134 using S3 data 
from GEO and the three LIGO interferometers. 
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